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The  transition  to  turbulence  in  a  boundary  layer  over  a  flat  plate  with  mild  surface 
undulations  is  simulated  using  the  parabolized  stability  equations  (PSE).  The  simulations 
incorporate  the  receptivity,  the  linear  growth,  and  the  nonlinear  interactions  leading  to 
breakdown.  The  nonlocalized  receptivity  couples  acoustic  perturbations  in  the  free-stream 
with  disturbances  generated  by  the  surface  undulations  to  activate  a  resonance  with  the 
natural  eigenmodes  of  the  boundary  layer.  The  nonlinear  simulations  display  the  influence 
of  the  receptivity  inputs  on  transition.  Results  show  the  transition  location  to  be  highly 
sensitive  to  the  amplitudes  of  both  the  acoustic  disturbance  and  the  surface  waviness. 
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I.  Introduction 


In  the  absence  of  significant  cross  flow  or  curvature,  transition  to  turbulence  generally 
results  from  the  amplification  of  traveling  eigenmodes  such  as  Tollmien-Schlichting  (TS) 
waves  or  Squires  modes.  These  traveling  modes  are  generated  through  the  process  called 
receptivity.  While  surface  vibrations  and  atmospheric  turbulence  can  activate  receptivity 
mechanisms,  experiments  have  shown  that  the  laminar-turbulent  transition  over  a  wing 
surface  is  strongly  affected  by  the  acoustic  field  generated  by  the  engines  and  the  turbulent 
boundary  layer  on  the  fuselage  [1]. 

The  interaction  of  the  acoustic  field  with  a  single  bump  has  been  the  focus  of  several 
investigations  in  the  past,  e.g.  Goldstein  h  Hultgren  [2]  and  Kerschen  [3].  These  analyses 
are  “decoupled”  from  the  downstream  disturbance  evolutions,  and  are  solved  independently 
to  provide  the  initial  amplitudes  for  the  evolution  calculations.  These  and  other  receptivity 
studies  have  been  primarily  theoretical  or  experimental,  but  there  are  some  ongoing  studies 
based  on  numerical  simulations,  such  as,  for  example,  Lin,  Reed  &  Saric’s  [4]  simulation  of 
the  receptivity  to  sound  for  a  flat  plate  with  a  blunt  leading  edge. 

A  new  mechanism  has  recently  been  forwarded  by  Crouch  [5,6]  based  on  small  nonlo- 
calized  irregularities  that  cover  the  entire  surface  (e.g.  undulations  or  roughness).  Possible 
origins  for  these  undulations  include  the  extrusion  and  plate  rolling  process  at  the  manufac¬ 
turing  stage,  and  heat  and  pressure  induced  distortions  during  flight.  Unlike  the  localized 
case,  the  receptivity  is  mixed  with  the  disturbance  evolutions  and  both  must  be  considered 
simultaneously.  Studies  indicate  that,  if  operative,  this  mechanism  is  likely  to  dominate  the 
receptivity  process  [7,8]. 

We  employ  the  parabolized  stability  equations  to  simulate  the  boundary-layer  transition 
over  a  wavy  surface;  we  include  both  the  receptivity  and  the  subsequent  linear  and  nonlinear 
evolutions.  The  relative  efficiency  of  the  PSE  simulations  accommodates  the  investigation  of 
the  effects  of  “receptivity  input  parameters”  on  the  path  to  transition.  These  investigations 
display  the  strong  effect  which  small,  i.e.  0(1),  changes  in  the  wall  roughness  height  and 
geometry  have  on  the  transition  process. 

II.  Wall  geometry  and  acoustic  field 

We  consider  a  flat  plate  immersed  in  an  incompressible  flow  field  with  acoustic  noi.se. 
The  plate  surface  is  covered  by  small  amplitude  irregularities  with  letigth  scab's  cotTiparabh' 
to  TS  waves. 

We  employ  a  Cartesian  coordinat<’  .sy.stem  with  the  average  plat('  .surfac<'  in  tin'  .r*-r* 
plane,  ;r*  measuring  streamwise  distance,  and  y"  the  distance  normal  to  tin'  i)late  (tin- 
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symbol  *  denotes  a  dimensional  quantity).  We  nondimensionalize  quantities  using  the  fn'e- 
stream  velocity  and  the  length  8*  =  '''*■*■* with  x’  being  the  location  where  the 

Reynolds  number  =  V^x*ju*  equals  one  million.  Results  are  |)resented  in  terms  of  the 
usual  stability  parameters,  R  =  \/R^  (distance)  and  F  =  10®(cJ*/27r)//*/(f ',*  )^  (freciuency). 

The  x-2-periodic  function  //(.r,r)  describes  the  surface  undulations.  We  represent  //  by 
its  Fourier  series  representation, 

'X>  oo 

f/(.r,r)=  ^  (1) 

r=— >x)  s=— 'X> 

where  the  coefficients  Wt,s  and  W_r,-s  are  complex  conjugates.  Additionally,  we  impose 
symmetry  in  thus  Wr,s  -  VVr,-.,.  We  use  values  of  |W|  of  the  order  of  0.002  (s('e  Table  1) 
which  translate  into  a  height  to  wavelength  ratio  of  order  1  /1 000.  At  STP  and  =  10  ?// /.s 
the  peak-to-peak  surface  variation  is  in  the  range  of  30  micromete’s,  while  at  =  100  m/s 
the  variation  is  about  10  micrometers. 

The  free-stream  acoustic  field  is  of  the  form 

OO 

u  =  .  V  =  lu  =  0.  (2) 

/=— OC' 

Associated  with  each  discrete  acoustic  freqtiency  At  is  a  velocity  field  having  a  Stokes  layer 
at  tlie  wall  satisfying  the  no  slip  boundary  condition,  and  matching  the  acoustic  field  in  the 
free-stream.  The  values  of  ui  in  (2),  and  of  cv  and  /3  in  { 1 )  represent  the  lowest  common  divisor 
of  the  set  of  frequencies  and  wave  numbers  present.  In  case  the  wall  spectrum  is  dominated 
by  sliari)  iieaks  but  the  acoustic  spectrum  is  not,  we  choose  u>  such  that  for  each  wall  mode 
{no,  kfi),  the  triplet  {lu;,  no,  k/i)  is  as  close  to  branch  I  of  the  TS  wave  neutral  stability  curve 
as  possil)le.  In  this  way  we  focus  our  attention  on  the  temporal-s|)atial  combination  that 
will  feed  the  greatest  amount  of  energy  into  the  eigenmodes.  Conversely,  if  a  flat  spatial 
si)ectrum  is  present  and  isolated  [)eaks  exist  in  the  acoustic  s})ectrum.  we  seh'ct  the  values 
of  fv  and  R  that  yield  triplets  close  to  branch  I. 

III.  Receptivity  mechanism 

Acoustic  disturbances  in  the  free-stream  generate  .Stokes  modes  within  llu'  lioundary 
layer.  In  th<’  incompressible  li(Tiit  th(>.se  modes  hav<‘onIy  temporal  ?7)odulation. 

V/ =  v/(,v)c''"' -f  c.r.  ,/  =  ...,-2.-1,  1.2 .  (3) 

.Vb'anwhile,  th('  truvui  flow  over  a  wavy  surface  produces  steady  wall  modes, 

=  +  r.c.  , 
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Figure*  ];  Rendition  of  surface  undulations  used  in  present  study.  The  normal  coordinate  is 
stretched,  (a)  (,'ase  “high”  and  “low”,  (h)  ('a.se  “rihlet”. 


=  ...,-2,-l.l,2,  ...  .  (4) 

These  modes  are  standing  waves  with  wave  numbers  given  directly  by  the*  surface.  Single- 
handedly,  neither  the  acoustic  nor  the  wall  velocity  fields  can  direct’y  energize  a  travel¬ 
ing  eigenmode  since  these  fields  lack  tlu*  necessary  spatial  or  temporal  variation,  respec¬ 
tively.  However,  the  simultaneotis  presence  of  both  fields  produces  traveling  waves  due  to 
the  quadratic  nonlinearity  of  the  Navier-Stokes  equations. 


V/.n./l- 


(■^) 


whose  form  is  identical  to  that  of  the  natural  eigenmode  of  the  l)oundary  layer,  except  the 
value  of  the  exponent  a  may  not  match  that  of  tlu-  eigenmode  of  the  boiindary  layc'r  stability 
equations,  which  has  the  form, 


^1,71, k  —  V; 


,l),  Jr., 


a,.  =  7,  +  ior- 


(b) 


The  receptivity  mechanism  i-^  illustrated  in  Figun*  2. 


The  key  ingredient  to  the  nonloralized  receptivity  process  is  a  resonance  which  results 
when  the  wall  wave  number  a  api)roaches  the  eigenmode  wave  number  Uf  —  7^  +  /Of  .  Here, 
7  is  the  growtli  rate.  Near  branch  I  the  growth  rate  is  small,  and,  for  an  appropriate  value 
of  o,  the  difference  |«f  —  oj  may  be  small.  The  resulting  respon.se  of  the  boundary  layer, 
under  the  forcing  provided  by  the  traveling  wave  may  be  quite  large.  The  detuning  \a,.  —  o| 


Acoustic  Wall 


Figure  2:  Components  of  the  nonlocalized  receptivity  model  (Spanwise  wavenumber  ftz  not 
shown). 

changes  as  the  modes  are  convected  downstream.  During  the  large  response  near  branch  I, 
energy  is  transferred  into  the  eigenmode  of  the  boundary  layer.  The  nonlocalized  receptivity 
model  [5,7]  shows  that  the  rate  of  energy  transfer  between  the  forced  wave  and  the  eigenmode 
is  proportional  to  the  rate  of  variation  of  the  forced  wave  response.  Farther  downstream, 
the  eigenmode  undergoes  the  typical  exponential  growth  characteristic  of  the  linear  regime. 
Receptivity  results  from  the  net  energy  transfer  into  the  eigenmode.  The  superposition  of 
the  forced  wave  and  the  eigenmode,  with  their  appropriate  wave  numbers,  provides  the  total 
(physical)  traveling- wave  disturbance. 

IV.  The  PSE  formulation 

As  a  consef|uenre  of  the  basic  flow  being  independent  of  the  spanwise  coordinate  r,  we 
can  ix'duce  tin'  nundu'r  of  unknowns  from  velocity  components  u,  v,  u;  (along  .r,  y,  and 
respectively)  and  i)ressure  p  to  oidy  u  and  v.  We  eliminate  pressure  by  taking  the  curl  of 
the  momentum  e(|uation,  a!id  ru  using  the  continuity  equation.  The  boundary  conditions  are 


zero  velocity  at  the  wall,  and  (except  for  the  acoustic  modes  and  the  mean-flow  distortion) 
vanishing  velocity  far  from  the  plate. 

The  parabolized  stability  equations,  commonly  abbreviated  to  PSE,  were  develojjed  by 
Herbert  &  Bertolotti  [9]  to  incorporate  nonlinearity  and  the  slow  growth  of  the  boundary 
layer  into  the  boundary-layer  stability  computations.  The  results  were  found  to  agree  with 
those  of  full  DNS  simulations  up  to  “’spike  stage”,  where  the  complexity  of  the  flow  rapidly 
spreads  beyond  the  resolution  of  the  PSE  code  [10,11,12]. 

Herein  we  e.xtend  the  PSE  forimilation  to  incorporate  the  nonlocalized  receptivity  model. 
We  express  the  velocity  field  in  a  series  in  time  (index  /),  in  x  (index  n)  and  in  c  (index  k). 

OC>  oo  (>0 

v(.r,(y,~,/)  =  Vb(j',.(/)  +  H  Z!  Z!  Vp(.c,?/.~,0  (7) 

/=  — oo  ?i=  — oo  A.—  -00 

where,  p  is  the  wave-vector  {l,n,k).  The  velocity  field  of  each  mode  is  partitioned  into 

Vp{x,y,z,t)  =  ylp(j;)vp(j:,y)xp(j:,c,f).  (8) 

The  function 

/Ip(ar)  =  exp  lp{s)ds  (9) 

incorporates  (in  a  sen.se  made  specific  below)  the  amplitude  of  each  mode  and  the  function 
X  incorporates  the  wavelike  part  of  the  mode’s  velocity  field 

\p(a:,  c,  f)  =  exp  [inaa:  +  ikf3z  —  (10) 

This  formulation  differs  from  the  regular  PSE  formulation  [12]  in  that  the  wavenumber  o 
in  (10)  is  not  a  function  of  ,r,  but,  rather,  is  held  constant  in  order  to  accommodate  the 
wall  Fourier  modes.  As  a  consequence,  the  maximum  allowable  step-size  yielding  acceptable 
accuracy  is  significantly  smaller  than  that  permitted  using  the  regular  forimilation.  The 
reason  for  this  reduction  is  that  differences  between  the  modes’  ])hysical  wavenumbers  (based 
on  u'^naxi  fhe  Fourier  wavenumbers  no  must  be  resolved  by  the  jirofile  functions 

Vp,  which  are  subject  to  the  “parabolization”  approximation  (12).  We  expand  on  this  topic 
later  on. 

The  acoustic  and  wall  modes  have  a  wave  vector  of  the  form,  respectively. 

p  =  (u;,0,0)  ,  (2u;,0,0)  ,  ...  ,  (Lu;,0,0) 
p  =  (0,o,0)  ,  (0,0, /^),  ...  ,  {Q,  No,  I\  ii). 

Both  the  forced  traveling  wave,  (5),  and  corresponding  eigenmodi'.  (0),  have  tlu'  same  fri'- 
quency  and  s|)anwise  wave  number,  and  approxiinatf  ly  the  sami'  streamwis*'  wave  numbc'r. 


Accordingly,  both  the  forced  and  the  eigenniode  can  be  resolved  by  the  Fourier  mode  having 
the  wave  vector 

p  =  {nu),noi,  kfj) 

lor  some  n  and  k.  In  particular,  the  eigenmode  wave  number  naf.{.r)  varies  with  x  and,  hence 
will  not  erpial  the  wall  wavenumber  no  except  possibly  at  one  point.  The  difference, 

Ayp  =  c‘”  (11) 

is  resolved  by  the  profile  function  Vp.  Thus,  the  dependence  of  Vp  on  x  in  eq  (4)  is  the  key 
factor  that  allows  the  PSE  to  capture  the  nonlocalized  recei)tivity  mechanism. 

The  partial  differential  equation  governing  the  velocity  profiles  Vp  is  obtained  by  substi¬ 
tuting  the  expansion  (3)  into  the  Navier-Stokes  equations.  For  the  streamwise  derivatives 
we  make  use  of  the  slow  change  with  .r  of  the  profiles  and  growth  rates  with  the  rule 


Up  Vp  +  mup  +  ~{m  ~  1  )ap 


P  ''tpi 


where  «p(.r)  =  7p(.r)  +  ina.  For  m  >  1  the  streamwise  derivatives  of  the  mean-flow  Vg 
are  zero,  in  accordance  with  the  boundary-layer  approximation.  Similarly,  in  (7)  we  drop 
second-  and  higher-order  derivatives  with  respect  to  x  of  Vp  and  Up.  Performing  harmonic 
i)alance  yields  an  infinite  set  of  coupled  partial  differential  equations  of  parabolic  type  in  x 
of  the  form 

Lvp  +  ~  12  Q  ’  ^P-r]  ’  (^•^) 

where  the  operators  L,  M,N  and  Q  contain  derivatives  with  respect  to  y  only.  The  sum¬ 
mation  on  the  r.h.s.  of  (3)  is  is  truncated  to  some  number  {L,N,l\)  in  the  numerical 
com|)utations.  Due  to  the  symmetry  in  ^  we  only  need  to  solve  for  modes  with  non-negative 
wave  numbers  in  /,  r,  and  c.  Upstream  traveling  modes  are  not  allowed. 

A  second,  and  closing,  set  of  equations  for  the  PSE  is  needed  to  define  the  values  of  7p 
as  a  function  of  r.  Thf'se  etpiations  remove  the  ambiguity  existing  in  (4)  wherein  both  Vp 
and  Ap  dep('nd  on  x.  In  the  present  work  we  use  an  integral  norm  (although  other  choices 
are  possibh') 

where  TZ  denotes  the  real  part,  and  f  denotes  the  complex  conjugate.  See  references  [10,1 1 ,12] 
for  more  details. 

We  transfer  the  z<'ro-slip  wall  boundary  condition  to  y  =  0  via  a  Taylor  series  exj)ansion 
about  y  =  0,  and  stop  the  expansion  at  terms  linear  in  \H\  since  this  quantity  is  assumed 


small,  although  including  higher  orders  in  \H\  can  be  done  in  a  straight  forward  fashion. 
Performing  harmonic  balance  yields  the  boundary  conditions  satisfied  by  Vp  for  each  p. 


Vp(;r,0)  = 


^VB(:r,0) 


+  E>^P- 


f^Vr(.C,  0) 


Ar{.r) 


The  initial  conditions  for  the  PSE  are  obtained  by  means  of  a  local  procedure  that  is 
composed  of  a  system  of  ordinary  differential  equations  and,  hence,  is  free  from  upstream 
influence,  as  described  in  [10,11,12].  Briefly,  the  local  procedure  is  derived  by  expanding 
Vp{x,y)  in  a  two-term  Taylor  series  about  the  x  location  of  interest,  and  collecting  terms 
of  order  0{\)  and  C?(A;c).  To  be  consistent  with  the  ordering  of  terms,  the  wavenumber  a 
and  frequency  u>  is  assumed  greater  than  order  0( /?“’).  C'onsequently,  modes  having  zero 
wavenumber  and  frequency,  such  as  vortices  and  the  mean-flow  distortion,  are  initialized 
differently,  as  described  below. 

The  xPSE  transition  analysis  tool-kit  has  been  employed  for  the  computations.  The 
partial  differential  equations,  (13),  are  transformed  into  algebraic  form  by  use  of  a  multi- 
domain  spectral  collocation  technique  in  t/,  and  a  finite  difference  discretization  in  x.  Five 
domains  are  used  in  y,  with  limits  at  [0,4],  [4,12],  [12,24],  [24,50],  [50,100]  and  u  and 
V  are  approximated  by  17  and  19  Chebychev  polynomials,  respectively,  in  each  domain. 
Asymptotic  boundary  conditions  are  imposed  at  y  =  100.  In  particular,  in  the  far-field  the 
velocity  fields  decay  exponentially  in  y  with  rate  ([  is  imposed,  using  (  =  (i  or  C2,  which  ever 
has  a  smaller  real  part. 


During  every  streamwise  step,  the  nonlinear  algebraic  system  is  solved  iteratively  by 
modifying  the  values  of  7p  until  the  normalization  conditions  (9)  are  satisfied  to  a  pre¬ 
determine  level  of  accuracy. 

A  further  element  affecting  accuracy  is  the  presence  of  (11)  in  Vp.  For  small  amounts  of 
separation  [oe  —  o|,  the  difference  in  wavenumber  is  well  captured  by  the  streamwise  change 
of  Vp.  For  larger  amounts,  however,  the  PSE  results  loose  accuracy.  A  more  detail  study 
of  this  issue,  as  well  as  the  effect  of  step-size  on  accuracy,  can  be  found  in  reference  [14].  It 
suffices  here  to  say  that  a  difference  in  [oe  —  cvl/cv  <0.1,  along  with  steps  sizes  Ax  <  5,  lead 
to  acceptably  small  errors.  Increasing  the  differences  above  10%  leads  to  an  underestimate 
of  the  growth  rates. 
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V.  Results 


The  presentation  of  the  results  is  divided  into  three  sections.  In  the  first  one  we  look  at  the 
effect  of  roughness  lieight  on  transition  using  the  surface  shown  in  Figure  la.  Then,  to  tins 
surface  we  add  a  spanwise  periodic  array  of  “riblets”  aligned  with  the  streainwise  direction,  as 
shown  in  fig  lb,  and  show  how  their  presence  triggers  transition.  The  initial  conditions  for  the 
steady  modes  in  the  riblet  case  need  special  treatment,  and  a  discussion  is  presented  prior  to 
the  results.  Afterwards,  we  investigate  the  effect  of  raising  the  acoustic  level  on  the  transition 
location.  Even  though  we  only  carry  a  few  modes  and  frequencies  in  our  calculations,  our 
results  lie  close  to  the  experimental  correlations  between  transition  location  and  free-stream 
turbulence  levels.  All  quantities  presented  herein  have  been  nondimensionalized  with  d*  and 
f/*  . 

OO 

To  specify  the  problem,  we  need  to  select  the  number  of  acoustic  modes  present,  and 
their  frequency.  Arbitrarily,  we  choose  the  frequency  F  =  56  and  its  subharmonic,  F  =  28. 
The  presence  of  the  subharmonic  is  needed  to  generate  the  traveling  wave  which  will  undergo 
parametric  resonance  with  the  fundamental  mode  at  some  point  downstream.  CJiven  these 
frequencies,  we  select  the  wavelength  of  the  wall  modes  to  maximize  receptivity.  We  choose 
the  fundamental  wavenumber  a  for  the  Fourier  decomposition  of  the  wall  undulations  to 
be  a  =  0.0871,  so  that  the  wall  mode  with  wavenumber  2a  will  closely  match  wavenumber 
Of,  =  0.1723  at  branch  I  for  a  TS  wave  of  frequency  F  =  56.  The  fundamental  spanwise 
wavenumber  3  is  .set  to  the  value  of  0.15,  so  that  3  ^  associated  K  and  H  type 

secondary  instabilities  are  nearly  maximally  amplified  [13]. 

While  our  selection  of  wavelengths  for  the  wall  modes  is  arbitrary,  we  note  that  within 
the  range  of  variations  we  are  dealing  with  (i.e.  10  to  30  micrometers)  most  wing  surfaces 
will  contain  a  large  selection  of  undulations  with  different  wavelengths,  amongst  which  our 
particular  choices  are  likely  to  be  present. 

The  spectrum  in  the  numerical  calculation  is  truncated  to  Z/  =  2,  N  =  2,  K  =  2  in  t,  ;r, 
and  except  for  the  riblet  case,  where  K  —  3  is  used.  The  integration  is  started  at  /?  =  300, 
which  is  sufficiently  upstream  to  include  any  significant  resonance  prior  to  branch  \  {R  =  589 
at  F  =  56),  and  marched  downstream  with  a  constant  step  size  of  A;r  =  3.5.  While  the 
triple  Fourier  sum  allows  for  27  distinct  modes,  many  of  these  modes  do  not  play  an  active 
role  in  the  transition  process  and  can  be  excluded  from  the  calculation.  For  example,  the 
modes  having  phase  speeds  of  either  2c  =  2u;/a  or  c/2  =  ujf2a  remain  at  roughly  constant 
amplitude  levels  throughout  the  computation,  and  their  presence  does  not  alter  the  results 
shown  herein.  Apparently,  the  phase  speed  difference  prevents  strong  participation  with  the 
interactions  jjropagating  at  the  speed  c.  Also,  our  results  show  that  the  interactions  between 
mofles  of  ef|ual  phase  speed,  but  different  from  c,  do  not  result  in  rapid  amplitude  growth. 


In  our  study  we  employ  three  different  plate  geometries,  which  we  label  “high”,  “low” 
and  “riblet”.  The  coefficients  in  (1)  are  real  and  have  the  following  values, 


Table  1 


high 

low 

riblet 

■1 

UjQliyfl 

QQQm 

BTTTirM 

1.0  xlO-3 

jsmsi 

0.0 

0.0 

1'iT‘ni‘M 

0.0 

0.0 

The  geometry  between  “high”  and  “low”  differs  only  in  the  peak-to-peak  variation  of  the 
undulations.  The  “riblet”  geometry  contains  two  additional  modes  which  describe  stream- 
wise  aligned  undulations,  similar  to  the  streamwise  “riblets”  used  in  turbulence  drag  reduc¬ 
tion,  only  that  our  “riblets”  are  not  sharp  peaked,  and  have  a  very  small  height. 


V.l.  Effect  of  roughness  height 

To  study  the  effect  of  surface  roughness  height  on  transition  we  consider  the  “high”  and 
“low”  geometries.  The  acoustic  modes  have  a  u  peak-to  peak  variation  of  0.0010  Voo^  (i-e. 
A\  =  Az  =  0.0005  in  equation  2)  while  the  peak-to-peak  variation  of  the  wall  modes  (Table  1) 
differ  amongst  each  other  by  at  most  a  factor  of  2,  hence  no  bias  exists  towards  one  particular 
mode.  This  scenario  contrasts  that  of  ribbon-induced  transition  where  a  two-dimensional 
mode  dominates  in  amplitude  prior  to  the  onset  of  secondary  instability. 

Figure  3  displays  the  amplitude  evolution  of  modes  in  the  “high”  roughness  case.  The 
transition  process  follows  the  well  known  subharmonic  route.  We  focus  our  attention  on  two 
modes;  the  two-dimensional  mode  (2,2,0),  which  develops  into  a  TS  wave,  and  the  (1,1,1) 
mode,  which  develops  into  the  subharmonic  mode.  Initially  each  mode  is  composed  solely 
of  a  forced  traveling  wave.  As  the  modes  propagate  downstream  the  nonlocalized  receptivity 
process  pumps  energy  into  the  eigenmodes.  Consequently,  the  (2,2,0)  and  (1,1,1)  modes  are 
of  nearly  the  same  amplitude  and  exhibit  similar  growth  rates  during  their  early  evolution. 
At  R  —  589,  the  (2,2,0)  mode  passes  through  branch  I.  The  eigensolution  component  of  the 
mode  undergoes  exponential  growth,  while  the  forced  traveling  wave  diminishes  downstream. 

The  development  of  the  subharmonic  mode,  (1,1,1),  can  be  analyzed  in  two  parts,  the  first 
one  at  locations  below,  say  R  =  1050  (including  the  sharp  dip  in  amplitude),  and  the  second 
one  for  locations  above  R  =  1050.  The  first  part  is  dominated  by  the  process  of  receptivity, 
while  the  second  part  is  dominated  by  the  parametric  resonance  with  the  (2,2,0)  mode.  The 
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F'igure  '.i:  F-foot  mean  square  of  maximum  ii.  velocity  versus  R  for  the  “higli"  roughness  case. 

(lip  in  amplitude  lU'ar  R  =  1000  is  caused  by  a  change'  in  phase  of  180  degrees  in  the  complex 
velocity  field  V(|,i  ij.  This  change  is  due  to  .r  dep('ndent  cliang('s  in  the  coefficients  of  the 
equations  whim'll  govern  the  iiiteractiou  betwe'en  tlu'  forced  mode  and  the  eigenmode.  At 
higher  levels  of  acoustic  forcing  the  receptivity  and  parametric  regions  overlap,  eliminating 
this  sharp  dip  in  amplitiuh'. 

The  (2,2,1)  mode  does  not  lead  to  fundanu'utal  (i.e.  K-tv})e)  dynamics.  The  “porpoising" 
in  amplitude  sc'eii  in  Figure  persists  evx'ii  with  a  four-fold  increase  in  the  acoustic  amplitude, 
wliich  increase’s  by  an  equal  amount  the  ami)litud('  of  the  traveling  modes,  but  only  slightly 
the  amplitude  of  the  steady  morh’s.  'flu’  ix’ason  for  this  lack  of  K-type  msonan'ce,  thus, 
cannot  be  exi)lained  simply  in  terms  of  a  thrc'shold  amplitude  of  the  (2,2,0)  mode.  Since  K- 
type  resonance  invedves  the  triad  interaction  between  the  (2,2,0),  (2,2,1)  and  (0,0,1)  modes, 
ruie  may  sus|)(’ct  an  unfavorable  phas(’  relation  in  the  triad  to  be  (pienching  of  the  resonance. 
A  four-fold  increase’  in  wall  mode’  amplitudes,  on  the  other  hand,  increases  the  amplitude’ 
of  all  mode’s,  including  the  (0,0,1)  mode,  and  the  flow  (lis])lays  a  mixed  ITtype  and  K-type 
transition. 

Additional  insight  into  the  eiu’rgy  transfer  between  traveling  modes  and  eigenmodes,  as 
well  as  phase  cancellation  giving  rive’  te)  the’  “por|)oising"  amplitudes,  e  an  be  obtained  from 
the’  perturbation  the’ory  [5]  wherein  the’  fe)re  e’el  trave’ling  wave  and  the  eige’iisolution  are’  kept 
elistinct.  In  the  PSF  formulatie)n  e)ne’  e  anne)t  e’asily  se’parate  the  se)lulion  inte)  the  sum  of  the’ 
twej  coiTipeenents. 

Lowering  the’  am|)lituele’ eef  the*  wall  unelidations  tee  the’  “hew"  h’ve’l  (se’e’  table’  1)  re’sults 
in  the  elynainies  sheewn  in  Figure’  T  Me’twe’e’ii  R  =  080  and  R  =  1200  the’  subharmonie' 


F  igure  4:  Root  mean  square  of  maximum  u  velocity  versus  R  for  the  “low"  roughness  case. 

undergoes  rapid  growth  as  in  fig  3,  however  the  amplitudes  of  the  (2,2,0)  and  (1,1,1)  modes 
remain  below  the  threshold  for  self-sustained  growth,  causing  these  modes  to  decay  after 
their  respective  branch  II  location  at  R  =  1116  and  R  =  1340.  The  (0,0,2)  mode  undergoes 
rapid  growth  past  R  =  1200  due  to  the  direct  forcing  from  the  ( 1,1 ,1 )-( 1,1 ,1 )  interaction. 
Further  downstream,  this  mode  persists,  albeit  slowly  decaying,  being  a  remnance  of  the 
earlier  parametric  resonance.  Similarly,  the  (0,0,1)  mode  grows  under  the  forcing  from  the 
(2,2,0)-(2,2,l )  interaction.  Botli  the  mean-flow-distortion  (0,0,0)  and  this  mode  undergo 
weak  algebraic  growth  past  R  =  1400  (see  section  V.2.1),  but,  eventually,  the  flow  returns 
to  its  undisturbed  state. 

The  results  of  the  “high"  and  “low”  roughness  cases  clearly  display  the  sensitivity  of 
the  transition  process  on  small  changes  of  wall  roughness  height  or,  equivalently,  on  small 
variations  of  acoustic  amplitude  levels.  This  sensitivity  highlights  the  difficulty  of  predicting 
transition  with  a  simple  criteria  such  as  c". 

V.2.  Effect  of  streamwise  undulations 

In  the  “high”  roughness  case,  displayed  in  Figure  3,  on(‘  can  ob.serve  the  birth  of  the 
mode  (1,1,0)  at  about  R  =  1170.  This  mode  is  of  interest  since  as  an  eigenmode  it  has 
an  unstable  region  that  spans  from  R  =  864  to  R  ^  2000.  This  mode  is  forced  via  the 
(0,0,1 )-( 1 ,1 ,1 )  mode  interaction,  and  since  the  wall  Fourier  mode  Wqij)  dirc’ctly  forces  the 
mode  (0,0,1)  through  the  wall  boundary  conditions,  w'e  looked  into  what  (dfi'ct  a  pt'iiodic 
array  of  streamwise  “riblets”  would  have  on  transition.  We  represent  tlu'se  “ribh'ts"  with 
two  Fourier  modes  with  wavehuigth  and  2/3.  (k)nse(|uently,  our  “ribh'ts"  are  more  blunt 
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and  smaller  than  those  used  in  turbulent  drag  reduction. 


V.2.1.  The  initial  condition 


The  generation  of  the  initial  conditions  for  the  (0,0,1)  and  (0,0,2)  modes  at  our  starting 
location,  R  =  300  posed  a  problem.  Several  runs  with  the  PSE  showed  that  the  evolution 
of  this  steady  mode  depends  strongly  on  the  initial  conditions  chosen.  For  some  conditions, 
algebraic  growth  ensued.  The  dependence  of  the  solution  on  the  initial  conditions  persists 
far  downstream,  hence  investigations  are  bia.sed  by  the  choice  of  initial  conditions  for  the 
steady  modes,  or,  equivalently,  on  the  birth  mechanism  of  these  modes.  Thus,  we  developed 
a  particular,  and  reproducible,  way  of  generating  the  initial  conditions.  Since  large  algebraic 
growth  can  obscure  the  dynamics  of  receptivity,  and  since  the  numerical  initial  conditions 
leading  to  such  large  growth  may  not  be  realizable  in  a  real  flow,  we  gauged  the  “correctness” 
of  our  initial  conditions  by  the  amount  of  algebraic  growth  present. 

We  generate  the  initial  conditions  by  extending  the  riblets  to  the  leading  edge  in  such 
a  way  that  a  separation  of  variables  approach  is  applicable  with  small  error.  The  riblet  is 
extended  linearly  from  its  full  height  at  the  PSE  starting  location,  x  —  xq,  back  to  zero 
height  at  the  leading  edge,  x  =  0, 

H(x,2)=  ^  0<x<x„ 

5=-2 


The  wall  forcing  (15)  involves  the  product  of  H  with  the  derivative  of  the  Blasius  U  compo¬ 
nent  of  velocity,  which  decreases  as  x“'/^,  therefore  the  forcing  (15)  increases  as  the  y/x  over 
0  <  X  <  xq.  Accordingly,  in  this  x  interval  we  introduce  the  following  form  for  the  (0,0,1) 
velocity  field, 

w(^,y,2)  =  u{y) 
v{x,y,z)  =  v{y)- 


X 


w{x,y,z)  =  w{y)—.  —e'  , 

Xq  V  Xq 


with  a  similar  form  for  the  (0,0,2)  field.  The  boundary  conditions  reduce  to  an  x  independent 
form, 

u(0)  =  -VWo.,r(0)  ,  u(0)  =  0  ,  th(0)  =  0 


where  /  is  the  .self-similar  stream  function  of  Blasius.  The  continuity  and  momentum  equa¬ 
tions  contain  coefficients  that  depend  on  x,  so  the  separation  of  variables  technique  is  not 
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strictly  valid.  However,  the  variation  with  x  is  slow  enough  to  be  negligible.  For  example, 
expanding  x  =  Xq  +  Ax,  the  continuity  equation  reduces  to, 

dv  .  Ax 

u  +  2- — h  2 — w  =  - - — u 

dy  xo  .ro  +  A.r 

The  right-hand-side  depends  on  x  but  is  0(R^')  near  x  =  .Cq.  VVe  set  the  right-hand-side 
to  zero  and  solve  the  equations  at  x  —  Xq  —  .300.  The  downstream  evolution  of  the  (0,0,1) 


Figure  5:  Amplitude  evolution  of  the  (0,0,1)  mode,  (a),  and  of  the  (0,0,0)  mode,  (b).  Dashed 
lines  show  cases  of  algebraic  growth. 

mode  is  shown  as  the  solid  line  associated  with  label  “(a)”  in  Figure  The  amplitude 
remains  nearly  constant  over  the  domain  of  integration,  showing  only  a  weak  transient  near 
/?o  =  300.  Negligible  algebraic  growth  is  present,  and,  in  this  respect,  our  initial  conditions 
successfully  generate  a  quasi-uniform  downstream  solution.  In  contrast,  the  dashed  line 
associated  with  label  “(a)”  corresponds  to  a  (0,0,1)  mode  generated  via  the  local  procedure, 
and  displays  the  linear  growth  (when  plotted  in  x  vs.  A)  characteristic  of  algebraic  growth 
[15,16]. 

The  mean-flow-distortion,  i.e.  mode  (0,0,0),  is  started  from  zero,  and  also  displays  al¬ 
gebraic  growth.  The  solid  line  associated  with  label  “(b)”  shows  the  development  of  the 
mean-flow-distortion  in  the  presence  of  the  (0,0,1)  mode.  The  growth  exhibited  is,  however, 
due  mostly  to  algebraic  growth.  This  fact  can  be  deduced  by  overlaying  on  the  Figure  the 
evolution  of  the  mean-flow-distortion  that  has  had  the  (0,0,1)  forcing  removed  after  a  few 
marching  steps  (shown  as  the  dashed  line). 
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V.2.2.  Effect  of  riblets  on  transition 

The  addition  of  riblets  to  the  ‘iow”  roughness  case  causes  the  flow  to  reach  the  later 
stages  of  transition.  The  amplitude  evohition  of  selected  modes  is  shown  in  Figure  (i.  In 
comparison  with  Figure  5,  the  (0,0,1)  mode  has  a  much  higher  amplitude,  and  the  (1,1,1)- 
(0,0,1)  forcing  generates  the  (1,1,0)  mode  at  about  R  =  1100.  While  initially  driven  mainly 
by  the  forcing,  the  eigenmode  component  of  (1,1,0)  undergoes  exponential  growth  (Branch  I 
is  at  /?  =  864)  and  becomes  the  dominant  component  past  R  —  1450.  Since  at  this  frequency 
the  eigenmode’s  total  amplification  AujAi  is  7208,  the  (1,1,0)  mode  grows  sufficiently  to 
produce  a  significant  back-forcing  on  the  (1,1,1)  mode  via  the  (1,1,0)-(0,0,1 )  interaction. 

Past  R  =  1650  the  (1,1,0),  (0,0,1)  and  (1,1,1)  modes  behave  like  in  the  vortex-wave 
interaction  study  presented  in  [10],  wherein  these  modes  eventually  lock  into  a  K-type  para¬ 
metric  resonance.  A  similar  triad  resonance  develops  between  the  ( 1,1,0),  (0,0,2),  and  (1,1,2) 
modes.  The  flow  at  /?  =  1970  reaches  a  strongly  nonlinear  stage  which  precedes  the  onset 
of  turbulence. 

We  like  to  refer  to  the  evolution  shown  in  Figure  6  as  the  “lateral”,  since  the  first  strong 
interaction  of  H-type  sets  the  stage  for  the  following  K-type.  This  progression  from  higher 
to  lower  frequencies  follows  the  slope  of  the  branch  I  of  the  neutral  stability  curve  as  the 
Reynolds  number  increases,  and  thus,  maintains  the  dominant  modes  near,  or  inside,  the 
eigenmode’s  region  of  instability. 


Figure  6:  Root  mean  square  of  maximum  u  velocity  versus  R  for  the  “Riblet”  roughness 
case. 


F'igure  7  gives  a  visual  synopsis  of  the  “riblet”  residts.  Shown  as  a  solid  line  and  a  dashed 
line,  respectively,  are  the  neutral  stability  curves  of  two-dimensional  and  three-dimensional 
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disturbances  given  b\'  the  linearized  stability  equations.  Superimposed  are  four  constant 
frequency  bands  that  represent  four  jjarticidar  modes  present  in  our  calculation.  The  bands 
display,  roughly,  the  location  and  extent  o(  the  receptivity  (light  gray),  linear  stage  (white) 
and  nonlinear  stage  (dark  gray)  in  our  computational  domain.  The  upper  of  tlie  two  bands 
at  F  =  of)  corresponds  to  th('  (2,2,0)  mode,  and  the  one  just  Ix'low  it  corresponds  to  the 
(2,2,1)  mod('.  Similarly,  tlu’  two  baiuls  at  F  =  28  correspond  to  the  (1,1,0)  mode  (above) 
and  the  (1,1,1)  mode  (below).  The  band  at  F  =  0  represents  the  (0,0,2)  mod(>.  For  this 
mode  the  receptivity  region  denot('s  tlu'  direct  forcing  by  the  VV(o,2)  vvall  mode, 

F=250 

2(J() 

l.-SO 

100 

50 


0 

0  500  1000  15{X)  R=2000 

Figure  7:  Visual  .synopsis  of  the  dynamics  shown  in  figure  6.  Integer  triplets  denote  the 
mode  represented  by  each  strip. 

V.3  Effect  of  acoustic  levels  on  Ruuhs 

We  close  this  section  with  the  results  of  our  investigation  into  the  effect  of  the  acoustic 
forcing  level  on  the  transition  Reynolds  number.  While  keeping  the  wall  geometry  fixed 
(“high”  case),  we  increa.sed  the  acoustic  amplitude  lewis  of  both  frequency  components 
F  =  56  and  F  =  28  from  At  =  0.0005  used  in  our  studies  above,  to  At  —  0.002,  0.01.  ami 
0.02.  The  corresponding  r.m.s.  amplitude  levels  are  \/2  times  larger.  Since  our  simulations 
do  not  reach  into  the  turbulent  flow  regime,  we  arbitrarily  defined  the  transition  Reynolds 
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iiiiinher  as  tlu"  location  wIhmc  tlic  sul)hannoiiic  motle  reaches  amplitude.  .At  this  point 
the  flow  is  engaged  in  the  strong  noidiiu'ar  interaction  charact(Tistic  of  spike  stag*',  and  liie 
turbulent  regime  usually  follows  within  a  coui)le  of  streamwise  wavelengths,  'lir/a.  W’e 
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F'igure  8;  Transition  Rj-  versus  free-stream  turbulence  level,  Tu.  Dark  bands  represent 
computetl  results. 

compare  our  results  with  the  summary  of  experimental  data  presented  by  Arnal  in  [1].  The 
experimeiitally  measured  (free-stream)  turbidence  level  was  based  on  the  r.m.s.  fluctuation 
of  all  three  velocity  components,  and  even  though  this  quantity  consisted  of  a  mixture  of  both 
sound  and  vorticity,  our  data  matches  qualitatively  with  the  experimental  ones,  as  shown  in 
Figure  8  (note  that  =  Rj.).  The  higher  numbers  of  Rxir  experimentally  measured  at  the 
cpiiet  limit  suggest  that  a  wall  smoother  than  the  one  modeled  here  was  used  in  the  wind- 
tunnel  tests,  while  the  lower  Rxtr  measured  at  higher  free-stream  turbulence  levels  suggest 
that  other  mechanisms  in  addition  to  those  considered  herein  are  at  work.  Considering  the 
simplicity  of  our  model,  the  “ball-park”  agreement  of  results  is  surprising. 
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VI.  Conclusion 


CJiveii  the  description  of  the  acoustic  field  spectrum  and  of  the  wavy  wall  geometry, 
we  have  followed  the  houudary  layer  response  through  the  regions  of  receptivity,  linear 
growth,  and  nonlinear  mode  interaction.  The  “‘bareness”  of  the  ini)ut  data  needed  for  the 
numerical  computation,  namely  the  W(„,*,.)  and  Ai  coefficients  in  e<|nations  (1)  and  (2). 
along  with  the  efficiency  inherent  to  the  marching  procedure,  facilitates  the  exploration 
into  the  effects  of  receptivity  parameters  on  transition.  Here  we  have  varied  (slightly)  the 
amplitude  of  the  surface  undulations  to  show  two  different  downstream  evolutions,  one  of 
which  reaches  spike  stage,  the  other  one  returning  to  the  Blasius  profile.  Adding  strt'amwise 
riblets  to  the  later  case  brings  about  spike  stage  through  a  ‘“lateral”  mechanism  involving 
the  transfer  of  resonance  from  higher  frerpiencies  to  lower  ones.  Our  results  on  the  effect  of 
acoustic  amplitude  level  on  the  transition  location  compare  qualitatively  with  wind  tunnel 
measurements  involving  free-stream  turbulence,  despite  the  simplicity  of  our  model.  Both 
the  acoustic  level  used  (0.007%  I’.m-s)  and  amplitudes  of  the  surface  waviness  (30/n/?. 

or  less)  make  our  investigation  relevant  to  flight  conditions. 
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